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We study fixed assemblies of touching spheres that can individually rotate. From any initial state, 
sliding friction drives an assembly toward a slip-free rotation state. For bipartite assemblies, which 
have only even loops, this state has at least four degrees of freedom. For exactly four degrees of 
freedom, we analytically predict the hnal state, which we prove to be independent of the strength of 
sliding friction, from an arbitrary initial one. With a tabletop experiment, we show how to impose 
any slip-free rotation state by only controlling two spheres, regardless of the total number. 


PACS numbers: 05.45.Xt, 89.75.-k, 45.40.Bb, 45.20.dc 

In a fixed assembly of touching spheres that can ro¬ 
tate individually, sliding friction at the contacts between 
spheres will generally slow down and finally stop the rota¬ 
tion of all spheres. However, if the assembly only contains 
even loops (bipartite assembly), it will instead drive to¬ 
ward a stationary slip-free rotation state which is free of 
sliding friction. This state happens to have at least four 
degrees of freedom. Here we introduce time invariant 
quantities and with them show, for assemblies with ex¬ 
actly four degrees of freedom, that no matter how many 
spheres belong to the assembly, one can always predict 
from the initial state of rotation the final state, which we 
prove to be independent of the type and strength of slid¬ 
ing friction. This allows imposing any slip-free state by 
only controlling two spheres, providing a method to con¬ 
trol the collective state of rotation of an assembly. With 
this work, we contribute to the understanding of the ro¬ 
tational dynamics of dense packings which are studied in 
the context of shear bands and seismic gaps [I1E1I5]- 
[9]. Since we propose a new way of controlling the rota¬ 
tion of spheres, mechanics and robotics is another area 
of applicability of our results. 

Consider three pairs of plastic spheres in contact al¬ 
ternatively displaced (at a stacking angle) on top of each 
other, with their positions fixed by ball rollers as shown 
in Fig.j^a). Forcing the two bottom spheres to rotate as 
indicated in Fig. [^b), we find a stationary state of rota¬ 
tion in which the top pair of spheres rotates more than 
three times faster than the bottom pair (video in the 
Supplemental Material). Since any two touching spheres 
have equal tangential velocities at their contact point, 
this state of rotation is slip-free and, therefore, free of 
sliding friction. This assembly is bipartite and as we 
show later, its slip-free state has four degrees of freedom, 
such that we call it a 4DOF assembly. Here, we first ex- 
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Figure 1. Operation of an assembly of spheres in a 
slip-free state, (a) Three pairs of plastic spheres in contact 
were stacked onto each other with an alternating stacking 
angle of 35 degrees and their positions were fixed using ball 
rollers. When forcing the bottom pair of spheres to rotate 
according to the scheme shown in (b) using external driving 
wheels, the other spheres adjust their rotation due to friction. 
In the stationary state the top pair of spheres rotates more 
than three times faster than the bottom pair. 


plore the kinetics of bipartite assemblies in general. We 
then show how to construct 4DOF assemblies, how to 
predict their final slip-free state from an arbitrary initial 
one, and how to control the rotation state experimen¬ 
tally. The discovery of the ability to control the rotation 
state of an assembly was unexpected. Our work is moti¬ 
vated by previous studies on bipartite assemblies of disks 
[IHZl fTOHT^ and spheres O |9] . 

Let us first consider two single spheres in contact that 
relax from an arbitrary initial rotation state toward a 
slip-free state, as shown in Fig. In the following, we 
color any sphere either red or yellow such that no spheres 
of the same color touch, what is only possible in bipar- 
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different tangential velocities equal tangential velocities 

Figure 2. Relaxation of a single pair of spheres toward 
the slip-free state. Rotating spheres with sliding forces 
at their contact (left) and in the hnal slip-free state (right), 
in which the contact moves along circles (dashed). Due to 
Newton’s third law of motion F 2 — — A, where Fi and F 2 
are the sliding forces acting on the hrst and second sphere, 
respectively. 

tite assemblies. We assume spheres are perfectly rigid 
and we only consider sliding friction. As long as the 
slip-free state is not reached yet, the tangential velocities 
of the spheres at the contact differ and the two sliding 
forces at the contact (one on each sphere) tend to re¬ 
duce this velocity difference. The two forces are oppo¬ 
site to each other and each force F produces a torque 
T = rxF, where r points from the center of the cor¬ 
responding sphere to the contact. The two torques are 
parallel and their magnitude is proportional to the radius 
of the corresponding sphere. Using the law of motion 
T = Id, we find independently of the type and strength 
of sliding friction 

h^i _ /-jN 

ri r2 

where U, I 2 , di, (^ 2 , ri, and r 2 are the moments of in¬ 
ertia, the angular accelerations, and the radii of the first 
and the second sphere, respectively. Note that since we 
assume that the centers of the spheres are fixed, the two 
sliding forces merely lead to torques which do not cancel 
each other such that angular momentum is not conserved. 

In bipartite assemblies with many spheres, a single 
sphere might have multiple contacts and multiple simul¬ 
taneously acting sliding forces. Each sliding force con¬ 
tributes to the angular acceleration of the sphere. We 
define df as the contribution to the angular accelera¬ 
tion di due to the sliding force at contact k, such that 
di = df, the sum running over all contacts. We write 
Eq. ^ analogously for a contact k between spheres i and 
j with the contributions d^ and d^. Eor simplicity, we 
consider the moment of inertia of a sphere i proportional 
to its mass rui and its radius squared, as, for example, 
for homogeneous solid spheres. We find 

miVid^ — mjVjd^ = 0 . ( 2 ) 

Summing this equation over all contacts k, we obtain 

'^Simiridi = 6, (3) 


where the sum runs over all spheres and is +1 if sphere 
i is red and —1 if it is yellow. Using Eq. (§ we define 

A :='^SimiriCoi, ( 4 ) 

i 

where ooi is the angular velocity of the sphere i. Equation 
^ shows that dAjdt = 0 , i.e., A is time invariant. 

We derive one further time invariant quantity. Let 
Xi be the position vector of the center of sphere i, 
where we choose the center of mass M of the entire as¬ 
sembly to be the origin of our coordinate system, i.e., 
M = 0. We multiply (dot product) Eq. ([^ by Xi and 
obtain miVid^ • Xi — rujrjdj • = 0. Since the vector 

Xij = Xj — Xi is perpendicular to d^, mjVjdj • Xij van¬ 
ishes. We subtract it from the previous equation and 
get miVid^ • Xi — mjVjdj • Xj =0. Analogously to the 

derivation of A, we thereby define the time invariant 
quantity B as 

B :='F SiVaiTiUJi • Xi- ( 5 ) 

i 

Later we use the quantities A G and R G R to 
predict the slip-free state of 4DOE assemblies. Now we 
identify how 4DOE assemblies can be constructed. The 
condition for the slip-free state of two spheres % and j is 
that their tangential velocities at their contact are equal. 
We formulate this as x ViXi^ = x r^x^i, where 
is the angular velocity of sphere % in the slip-free state 
(sf), and Xi^ is the unit vector pointing from sphere % to 
sphere j. We rewrite this condition in analogy to Ref. [8] 
as 

SjrjLof - SiTiLdf = CijXij, (6) 

which defines c,j S R that uniquely relates iof to 
The slip-free state is uniquely defined by the angular ve¬ 
locity of a single sphere i and the parameters Cij of all 
contacting spheres i and j. Since ujf G R^, the number 
of DOEs of the slip-free state is equal to three plus the 
number of DOEs of the set of q^’s, such that we only 
have a 4DOE assembly, if the set of q^’s is restricted to 
a single DOE. Therefore a single pair of touching spheres 
with a single Cij is the simplest 4DOE assembly. Eor any 
longer open chain of spheres, we have an independent 
Cij for each contact as shown in Eig. I^a), not resulting 
in a 4DOE assembly. Eor a bipartite loop with an even 
number of spheres A" > 4, the q^’s of the contacts are 
not independent and have N — 2 DOEs in case all the 
centers of the spheres are coplanar, and N — 3 otherwise 
(details in appendix |^. Therefore, the only loop that is 
a 4DOE assembly is a loop of four spheres whose centers 
are not coplanar (see e.g. Eig. [^b)), which we denote as 
a noncoplanar-4 loop, where all q^’s have to be equal 
to a single parameter c. Starting with a pair of spheres, 
every other 4DOE assembly can be constructed by it¬ 
erative extension in the following ways, as long as the 
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Figure 3. Parameters describing the slip-free state. 

The Cij parameters that uniquely relate the angular veloci¬ 
ties of contacting spheres i and j in the slip-free state are 
independent for an open chain (a) and are all equal to one 
single c for a noncoplanar-4 loop (b), which is a 4DOF as¬ 
sembly. For Cij = 0, the angular velocities of spheres i and j 
are antiparallel. 


assembly is bipartite. Either connect two 4DOF assem¬ 
blies by establishing two contacts involving two spheres 
of each assembly such that the centers of the four spheres 
involved are not coplanar. Or connect a single sphere via 
two contacts to a 4DOF assembly such that the centers 
of the involved spheres are not collinear. With these op¬ 
tions as minimal requirements, one can even form the 
space-filling assembly presented in Ref. [8] . Both ways of 
extension ensure that the q^’s of all contacts are equal 
to a single c (details in the appendix [B] and [C|) , such that 
given Eq. the slip-free states of 4DOE assemblies can 

be uniquely described by 

wf = ^ + cxA , (7) 

where G is a vectorial reference quantity. 

We can now predict the final state toward a 4DOE 
assembly will drive by using Eqs. Q, ([^, and Q (details 
in the appendix to find 


UJ. 


sf 


_ Sj / A B ^ 

vAm^ 


(8) 


with M = and I = Equation ^ 

shows that the final slip-free state can be predicted using 
the time invariant quantities A and B defined in Eqs. 

and <§• Surprisingly, the final state is independent 
of the type and strength of sliding friction. Different 
sliding forces merely lead to different kinetic pathways 
toward the slip-free state as shown in the appendix |E| 
Eigure pictures the role of A and B regarding the 
slip-free state predicted by Eq. (§. Eor R = 0 all axes of 
rotation are parallel to A. Eor B 0 they all intersect 
at X = IAj (BM), and a sphere located at x would be at 
rest and would rotate faster the larger its distance to x 


Remarkably, if one blocks one sphere, generally not all 
spheres will stop rotating but the assembly will instead 
relax toward a slip-free state where all axes of rotation 
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Figure 4. Description of different slip-free states. 

4DOF assembly in different slip-free states. For A 0 and 
R = 0 all axes of rotation are parallel to A. For R 7^ 0 all axes 
of rotation of the spheres meet at the position x = IA/(BM), 
and in particular for A = 0 they meet at the center of mass. 


intersect at the center of the blocked sphere (details in 
the appendix Blocking two or more spheres at the 
same time will force every sphere to stop. Assemblies 
with more than four DOFs allow more than one sphere 
to be blocked without stopping all spheres. 

To impose any desired slip-free state of a 4DOF as¬ 
sembly within our model, one can apply instantaneous 
changes to two arbitrarily chosen angular velocities as 
we show in detail in the appendix 0 where we addi¬ 
tionally show how to determine the total mass M, the 
center of mass, and the parameter I by accessing only 
two spheres (appendix 0- A real system that is sub¬ 
ject to e.g., rolling, torsion, and air friction would always 
come to rest from any initial state. Nevertheless, one can 
preserve any desired slip-free state of Eq. (compare 
Fig.§ by preserving two angular velocities accordingly. 
Ffow close the real stationary state will reach the desired 
state depends on various details of the assembly such 
as spatial arrangement, material, size, contact forces be¬ 
tween spheres, and the fixing structure. We conducted 
a simple experiment as a first demonstration on how the 
reached stationary state might depend on the spatial ar¬ 
rangement. Two horizontally contacting pairs of hollow 
TPU spheres were stacked onto each other as shown in 
Fig. B a) with varying stacking angle a (compare Eig. 
I^c)), forming a noncoplanar-4 loop (4DOE assembly) 
for 0° < Of < 90°, and a coplanar-4 loop (5DOE assem¬ 
bly) for a = 0°. We fixed the positions of the spheres 
using ball rollers, allowing the spheres to rotate. The 
two lower spheres were forced to rotate with equal ab¬ 
solute angular velocity around the axis connecting their 
centers. When we force them to rotate in opposite di¬ 
rections as shown in Eig. [^c), we impose a state where 
all rotation axes are parallel and the absolute angular 
velocities are equal, such that we call this the constant 
mode. When forcing both to rotate in the same direc¬ 
tion as shown in Eig. |^d), we impose for a > 0° a state 




4 



2.0 


amplifying mode 


1.Q f 


constant mode 


0 20 40 60 80 

a (degrees) 




amplifying mode 


Figure 5. Tabletop experiment to compare with 

the model, (a) Two pairs of contacting TPU spheres were 
stacked onto each other and their positions were fixed using 
ball rollers. The two lower spheres were forced to rotate with 
equal absolute angular velocity using driving wheels attached 
to motors, (b) Ratio of the angular velocities of an upper and 
a lower sphere (experiment (symbols) in comparison with the¬ 
ory (lines)) versus different stacking angles a (compare (c)), 
for two different modes of rotation indicated in (c) {constant 
mode) and (d) {amplifying mode). 


where all axes of rotation meet at the contact between 
the two lower spheres. Then the two upper spheres will 
rotate at a higher angular velocity as we have shown in 
Fig.[T] therefore, we call it the amplifying mode. Figure 
l^b) shows the measured ratio of the angular velocities 
of an upper and a lower sphere (symbols) as well as the 
prediction by our model (lines). In the constant mode the 
experiment agrees well with our prediction for all angles 
a. In the amplifying mode, we find excellent agreement 
for a « {60°, 84°}. For a « {12°, 24°, 36°}, the angular 


velocities of the upper spheres were on average smaller 
than the predicted ones, and we did not observe smooth 
but partially jerky rotations becoming jerkier for decreas¬ 
ing a. In the amplifying mode, the spatial arrangement 
for low a triggers this deviation which originates from 
the elasticity of the spheres. In our experimental setup, 
elastic deformation can lead to jerky motions by causing 
displacements of spheres and by the formation of non- 
negligible contact areas between spheres that can lead to 
torsional forces as we explain in detail in the appendix 
|Tj Experimental details can be found in the appendix 
and videos in the Supplemental Material. 

This Letter advances the understanding of the rota¬ 
tional degrees of freedom of assemblies and thus of the 
internal dynamics of shear deformation and seismic gaps. 
In a technological perspective, the ability to control the 
rotation state of an assembly of rotating spheres in con¬ 
tact is a newly discovered functionality with a general 
but yet unexplored potential. It is likely to find use in 
mechanics and robotics to control the orientation and 
rotation of spheres. The possibility to amplify the angu¬ 
lar velocities of spheres along a certain contact network 
could be employed as an alternative to power transmis¬ 
sion gears. Since the ability to control the rotation state 
is robust against changes in the spatial arrangement, as 
long as contacts are conserved, assemblies allow for de¬ 
sired or inevitable displacements during operation. 
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Appendix A: Degrees of freedom of c^^’s for loops 


Let us consider a bipartite loop with an even number 
of spheres N > 4. The angular velocities and oof of 
two contacting spheres i and j in the slip-free state (sf) 
are uniquely related to each other via a parameter Cij 
according to 

SjrjLof - SiViLdf = CijXij, (Al) 


where is the radius of the sphere i, is +1 if sphere i is 
red and —1 if it is yellow, and Xij = Xj — Xi is the vector 
pointing from the center Xi of sphere i to the center Xj of 
sphere j. Since for every contact Eq. (Al) needs to hold, 
one finds the constraint 


C12X12 + C23X23 + . . . + CniXnI = 0 , 


(A2) 


where the indexes 1 to A are given to the spheres in con¬ 
secutive order. The set of q^’s has N — R DOFs, where 
R is the rank of the 3 x A matrix {xi 2 X 23 • • • ^ is 

equal to the number of linear independent columns of the 
matrix and is either two, when all centers of the spheres 
are coplanar, or three ot herw ise. For noncoplanar-4 
loops, it follows from Eq. (A2) that all q^’s need to be 
equal and can be replaced by a single parameter c. 


Appendix B: Ensure that all c^j’s are equal 


To extend a 4DOF assembly, one needs to ensure that 
all Qj’s are equal to a single c. Starting from any 4DOF 
assembly, one finds for any pair of contacting spheres i 
and j that 


SjrjLof - SiTiLuf = CXij. 


(Bl) 


Let us first relate the angular velocities of two spheres 
i and k which are not in contact, but are both touching 
a third sphere j. We sum Eq. (Bl) and an analogous 


equation for the contact between sphere j and k to find 

SjrjUjf - Siniof + SkrkCol‘' - SjrjLof = cXij + cxjk, 

- SiViCof = CXik, 

(B2) 

which shows the same relation as Eq. (Bl) for the non¬ 


contacting spheres i and k. Therefore, Eq. (Bl) is valid 
for any pair of spheres in the assembly. 

There are two options of how to extend a 4DOF as¬ 
sembly. One way is to connect two 4DOF assemblies A 
and B. Before they are connected, each of them has an 
independent parameter c which we call ca a nd that 
describe the slip-free state according to Eq. ( |B1[ ), which 
holds for any two spheres i and j in the assembly. If we 
want to form another 4DOF assembly by connecting A 
and B, we need to make sure that the way of connecting 
enforces ca = cb for the slip-free state. This can be done 
by involving two spheres of each assembly to couple the 



Figure 6. Ways of extending 4DOF assemblies. Ev¬ 
ery 4DOF assembly can be iteratively constructed using two 
rules of connecting as minimal requirements. Either connect 
two 4DOF assemblies A and B by establishing two contacts 
involving two spheres of each assembly such that the centers 
of the four spheres i, j, k, and I involved are not coplanar 
(a). Or integrate a single sphere /c in a 4DOF assembly A by 
establishing two contacts between them such that the centers 
of involved spheres are not collinear (b). 


two parameters ca and c^- To involve two spheres of 
each assembly we need to establish at least two contacts. 
If we connect sphere i and j from assembly A to sphere 
k and I from assembly B as shown in Fig.a), and con¬ 
sider the fact that A and B are 4DOF assemblies with 
parameters ca and c^, respectively, we can establish the 
constraints 


- SiTiLof = CAXij, 

(B3) 

SinCof - SjTjUjf = CjiXji, 

(B4) 

Skfki^k - smojf = CBXik, 

(B5) 

SiViuf - SkVkOJk = CkiXki, 

(B6) 

which if combined lead to a single constraint 

^A^ij T CjiXji “h CbXi]^ “h — ( 

). (B7) 

Only if one finds three linear independent vectors among 
Xij^ Xji^ xik^ and Xki^ i.e., if the centers of the four 
spheres i, j, /c, and 1 are not coplanar, Eq. (B7) enforces 
Ca = Cji = Cb = Cki, such that the resulting assembly is 
a 4DOF assembly. Note that spheres i and j of the bi¬ 
partite assembly A can be of any color, as long as the 
final assembly is also bipartite. Another way to extend 
a 4DOF assembly is to integrate a single sphere k by es¬ 
tablishing contacts to two spheres i and j of a assembly 
A, as shown in Fig. |^b). The constraints 

SjVjUof — SiViUof = CAXij ^ 

(B8) 


(B9) 

SiTiuf - SkVki^of = CikXki, 

(BIO) 

can be added to obtain 


^AXij T CjkXjk T C-ikXki — 0. 

(Bll) 
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Only if the c enter s of the spheres j, and k are not 
collinear, Eq. (Bll) enforces ca = Cjk = Cik- In that case, 
the resulting assembly is a 4DOF assembly. Starting from 
two single spheres being the most simple 4DOF assembly, 
one can construct every other 4DOF assembly iteratively 
using the two presented rules for connecting as minimal 
requirements. The formation of additional contacts dur¬ 
ing the process does not change the fact that the resulting 
assembly is a 4DOF assembly as long as it is bipartite. 


Appendix C: Description of the slip-free state for 
4DOF assemblies 



. magnitudes of 
sliding forces: ^q“ 2 
-prop, to -4 

Vj-Vi 

-constant 10 ^ 



-h- 

time 3 


slip-free state 

I 

I - > 


Since Eq. (Bl) holds for any pair of spheres i and j 
in the assembly, we can express the angular velocity of 
any sphere j as a function of the angular velocity of a 
reference sphere i as 


= — (sinuf + cXij). (Cl) 

To describe the slip-free state in a general way, we choose 
the reference sphere i not to be an actual sphere of the 
assembly, but an imaginary reference sphere with angular 
velocity ujf = f], = 1, and = 0 to obtain 


Figure 7. Different sliding forces merely lead to 
different kinetic pathways toward the slip-free state. 

4DOF assembly that relaxes from randomly chosen initial an¬ 
gular velocities toward a slip-free state with sliding forces 
with magnitudes proportional to the difference in tangen¬ 
tial velocities (solid) and with constant magnitudes (dashed). 
A = \ui — a;f^| measures the deviation from the slip-free 
state. Inset: Exponential decay of A found for the propor¬ 
tional force. The final slip-free state is independent of the 
sliding forces in contrast to intermediate states with same A. 
Simulation time in seconds with a time step of s. 





(C2) 


time invariant quantities only as 


Appendix D: Predicting the slip-free state using A 
and B 



(D3) 


The slip-free state as written in Eq. (C2) can be pre¬ 
dicted using the time invariant quantities A and B de¬ 
fined as 


i 


SifniriCui 


(Dl) 


Appendix E: Different sliding forces lead to different 
kinetic pathways 


and 


i 


SiTriiriUJi • Xi. 


(D2) 


We write Eq. (C2) for a sphere i, multiply it by Sifriiri, 
and sum over all spheres to end up with A on the left 
hand side and to eliminate c from the right hand side, 
because J]) • crriiXi = cM = 0, since we defined the center 


of mass M of the entire assembly to be the origin of 
our coordinate system (M = 0). We find ^4 = A/M, 
with the total mas s M = Second, we multiply 

(dot product) Eq. ( |C2[ ) written for sphere i by SimiViXi 
and sum over all spheres to end up with B on the left 
hand side and eliminate Q on the right hand side to find 
c = H/1, where I = mi\xi\‘^. We then formulate the 

angular velocities of the slip-free state as a function of 


The final slip-free state is independent of the type and 
strength of sliding friction. Different sliding forces only 
lead to different kinetic pathways toward the final state 
as illustrated in Fig. which in general also depend on 
the geometry of the assembly and the moments of in¬ 
ertia of the spheres. We show a noncoplanar-4 loop, a 
simple 4DOF assembly, relaxing from a random initial 
confi guration toward the slip-free state predicted by Eq. 
(D3) for two different types of sliding forces. Any sliding 
force Fij acting on sphere j due to contact with sphere 
i points in the direction opposite to the relative velocity 
Vij = Vj — Vi at the contact. We considered two cases. 
First, Fij = —avij, where the magnitude of each force 
is proportional to the difference in tangential velocities, 
with <j = 3 and second, Fij = —avij^ where the magni¬ 
tude is constant, with a = 0.05 and with Vij being the 
unit vector along Vij. 
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Appendix F: How blocking a single sphere influences 
the slip-free state 


From Eqs. ( |D1[ ), ( |D2| ), and (D3), one can derive the 
new A (Anew) and B (Bnew) from the previous (old) val¬ 
ues (Aoid and Bold) when blocking sphere i as described 
in the following. By applying a perturbation Auji to the 
angular velocity of a sphere i one can impose the changes 




SiTTlir, 


■Auji 


and 


AiB = 


SiTriirilXuJi 


iAi^i • Xi. 


(Fl) 


(F2) 


Remember that the origin of the position vector Xi of 
the center of any sphere i is the center of mass of the 
entire assembly. Blocking a single sphere permanently 
and letting the assembly relax to a slip-free state has the 
same effect on A and B as applying a perturbation Aooi 
that leads to the sphere i being at rest in the slip-free 
state, i.e., = 0. To ensure Acl;®^ = 0 we find from 


Lof = ^ 


M 


+ 


Br, 


= 0 , 


that 


^new -MBnewXi/L 


(F3) 


(F4) 


Fro m th e fact that AiA • Xi = AiB (compare Eqs. (Fl) 
and (F2)), we find using A^A = Anew — ^oid and AiB = 
old') that 


^new 


^new — 


Bold — A 


old ’ 


l + M\xX/I 


(F5) 


Appendix G: Imposing a desired slip-free state 


We want to impose a desired change AA and AB by 
applying external changes Auji and Aujj to the angu- 


to the center of sphere j and subtract Eq. (G2) from it 
to find 


cj,: = 


AA • Xj — AB 

SjmjrjAujj • Xj 


(G3) 


such that we ca n choose any AuJi as long as AuJi-Xij ^ 0. 
From Eq. (G3) we know Auji since Auji = uJiAuji and we 
obtain using Eq. (|^ that 


Aujj = (AA — . 


iA(Xi)/{6 


,)• 


(G4) 


Appendix H: How to determine the total mass M, 
the center of mass, and the parameter I by accessing 
only two spheres 

It is possible to determine all global quantities rele¬ 
vant for the final slip-free state by accessing not more 
than two spheres o f a 4 DOF assembly regardless its size. 
Compared to Eq. (D3), relevant are the total mass M, 


the center of mass, i.e., the origin of the position vec¬ 
tors Xi, and the parameter I. Starting from any slip- 
free state, one can separately apply two changes Aiu;^ 
(first) and A 2 (Ai (second) to an accessible sphere i. One 
needs to wait after each change till the slip-free state is 
reached, and determine the corresponding changes Aiu;®^ 
and A 20 jf in angular velocity from the previous to the 
new slip-free state. Applying a change Auji to a sphere i 
leads to a change Ad;®^ of the angular velocity of sphere 
i between the previous to the new slip-free state. We use 


A ^sf ( A^A AiB ^ 


(HI) 


and Eqs. and ( |F2[ ) to find 

Acjf = mi(ALOijM -h (Acji • Xi)xill). (H2) 


Let us derive two scalar equations from Eq. ( |H2| ) by first 
squaring it (dot product) to find 


lar v^ity of sphere i and j, respectively. Eqs. ^ \AOf\‘^/m‘f = \AOi\^/M‘^ ^ 2{AOi • Xi)^/(MI)^ 


and (F2) show the changes A^A and AiB imposed by 
the change Auji. The individually imposed changes in 
A and B need to sum up to the desired change. 


(H3) 


I.e., 


A.A and AB = AiB 


A,B. 


AA = AiA + 

Auji = uJiAuJi _ _ _ 

change Auji. From Eqs. (|F1[) and (F2) we find 


We define 


where uji is the unit vector of the external 


A A = AiA -f AjA = SirriiriUJiACji + SjrrijrjAujj (Gl) 
and 

AB = AiB AjB = SiTriiriUJiACji • Xi + SjirijrjAujj • Xj. 

(G2) 


We multiply (dot product) Eq. (Gl) with the position 
vector Xj pointing from the center of mass of the assembly 


(AtJj • Xi)‘^\xi\‘^/I^ 

and second multiplying it (dot product) with Auji to find 

Acjf • Aui/rrii = \Aui\^/M + (AcJi • Xi)‘^/L (H4) 

Combining the two we can eliminate the term (AOi • Xi)^ 
and after some rearrangements we find 

1/M2 + (2/M + \xi\^/I)iACofACoi/imi\ACoi\^) - 1/M). 

JH5) 

For a singl e ch ange Auji and its induced change Ad;®^ we 
see in Eq. (H5) that we have the two unknown quantities 



















M and \xi\^/I. So with two changes AiWj and A 2 Wi and 
their induced changes Aia}|^ and A 2 w|^ one can find 


M = {b+ y/b^ - 4ac) /(2a), 
a = \f2g1 - /152I, & = 1/2 - fi\, c = \g2-gi\, 
U = |A„cJf|V(mf|A„cDi|2), 
gn = A„wf ■ A„Wi/(mi|A„Wi|). 


(H6) 


To locate the center of mass of the assembly we can 
use Eqs. (iDll), (|D2|), and (|D3|) to find a vector 


— rriiAuji/M^ 


(H7) 


that is parallel to the vector Xi pointing from the center of 
mass to the center of sphere i. By applying an additional 
chan ge A ooj to a second accessible sphere one can use 


Eq. (H7) to obtain a vector x', parallel to Xj . The center 


of mass is the only point in space that can be reached 
both going along x[ from the center of sphere i and going 
along x'j from the center of sphere j, in the general case 
wher e x[ is not parallel to x'j. At last, we can solve Eq. 
I4|) for / to find 


I = {Auji • Xi) / (Acjf • Auilrrii - \ALUi\‘^/M) . (H8) 


Appendix I: How elasticity can lead to torsional 
forces 

As in the amplifying mode of operation shown in Eig. 
ia), elastic spheres can lead to torsional forces due to 
the formation of finite contact areas between spheres as 
shown in Eig. [^b). Let us assume that the velocity of 
each contact point is proportional to its distance to the 
rotation axis of the sphere. If not all points of the contact 
area have the same distance to the rotation axis, the 
contact area shows an inhomogeneous velocity profile. 
This is the case for the contact areas between lower and 
upper spheres in the amplifying mode as shown in Eig. 
[^c). The angular velocities are assumed to be the ones 
predicted by our theory and in this example, we have 
a = 36°. To get a better view on the inhomogeneity of 
the velocity profiles, we plot the difference to the velocity 
in the center of the contact areas as shown in Eig. I^d). 
There one can see that the stronger inhomogeneity of the 
upper contact area will lead to forces on both contact 
areas as shown in Eig. [^e), that sum up to a torsional 
net force acting on the spheres. 


Appendix J: Experimental details 

We used hollow TPU spheres with the trade name Su¬ 
per High Bounee Jump Ball with a 10 cm diameter. The 
positions of the spheres were fixed using bolt fixing ball 
transfer units (ball rollers) with 19 mm carbon steel balls 
manufactured by Alwayse Engineering Limited. Eor each 
sphere two to three ball rollers were placed such that with 



Figure 8. Elasticity can lead to torsional forces, (a) 

Amplifying mode of operation: upper spheres (3,4) rotate 
faster than lower spheres (1,2). (b) the elasticity of spheres 
leads to a contact area (striped) between spheres instead of 
point-like contacts assumed in our theory, (c) Velocity pro¬ 
files of the contact areas between sphere 4 (top) and sphere 
1 (bottom), (d) Profile of the velocities after subtracting the 
velocity at the center of the contact areas, (e) Resulting forces 
acting on the contact areas lead to a torsional net force on 
the spheres (length of arrows are proportional to the velocity 
difference between touching points). 


all spheres inserted, all positions stay locked in operation. 
On the thread of the ball rollers, we placed springs with 
37.4 N spring rate followed by a short brass pipe and a 
screw nut, such that the ball rollers could be fixed at the 
pipe, using boss head clamps and rods, and their posi¬ 
tions could be slightly adjusted using the screw nut. In 
the case of too large contact force between the ball roller 
and a sphere, the spring gives in. Due to the usual appli¬ 
cation of ball rollers in cargo transfer, they come greased. 
Since we want to prevent the spheres from slipping on 
another, we degreased the ball rollers using petroleum. 
They were one by one dipped into a cup of petroleum 
such that the ball compartment is covered. Each ball 
roller was gently patted and rolled on the bottom of the 
cup for about one minute. Then, before they were dried 
on air, the petroleum was shaken off and absorbed with 
household paper. All ball rollers were cleaned in this 
way four times, where the petroleum was replaced after 
having washed every ball roller once in it. 

To control the rotation of the two lower spheres, we 
used for each of them a 31.7 mm large foam rubber tire 
that we fixed, using a propeller hub, on a gear motor with 
nine turns per minute and a peak torque of 1.8 N when 
operating at 12 V. The positions of the wheels were fixed 
such that the contact forces between the wheels and the 
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spheres was large enough to prevent slip during operation 
and small enough to prevent strong deformation of the 
spheres. 

Every operating setup was caught on video over five 
minutes. The angular velocities were measured by man¬ 
ual video analysis. Signs were drawn onto the spheres in 
advance to ease video analysis. We measured the angular 
velocity as the number of times a sphere turned during 
the five minutes divided by the exact number of passed 
seconds. We assumed the rotation axes of the spheres to 
be fixed over time, which turned out to be not always 


exactly the case especially for low values of the stacking 
angle a, but a reasonable assumption. 

To stack two or three pairs of spheres for a rather sim¬ 
ple demonstration, the hollow TPU spheres and the ball 
rollers fixed in springs worked well. To build assemblies 
with more spheres, the TPU spheres will turn out to be 
too soft at some point, and due to compression in dif¬ 
ferent directions depending on the operation mode, they 
will not suit anymore. For optimal slip-free operation, we 
suggest spheres to be as light and stiff as possible while 
having as much grip between spheres as possible. 



